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We prove, by realistic microscopic calculations within the sp'^s* Tight Binding method for GaAs 
(110) and (100), that the surface optical properties are not influenced by long-range crystal termi- 
nation effects, and hence that they can be consistently studied considering slabs of limited thickness 

... ( 20 - 30 j4). The origin of derivative-like and bulk-like lineshapes in Reflection Anisotropy Spectra 

0^ \ is also discussed, analyzing the effects arising from possible surface-induced reduction, broadening, 

Oi ■ and shifting of the bulk spectrum near the surface. 
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(T^ ! I. INTRODUCTION 

^. 

Reflectance Anisotropy Spectroscopy (RAS) and Sur- 
face Differential Reflectance (SDR) are surface-sensitive 

vQ : optical techniques, and are used to obtain information on 

Tjj- the atomic and electronic structures of surfaces j|]. In 

Cn . the early times of these spectroscopies, the spectra were 

^^ ' explained in terms of transitions across surface states, 

^^ ] although this view has been contrasted by calculations, 
showing that surface-geometry effects could also deter- 
mine the spectra through surface perturbations on the 

'S optical matrix elements of transitions across bulk states 
r^ 0. Now the attitude seems to be reversed: after real- 

'T' : izing that many RA lineshapes are similar to the imagi- 

'"^ nary part of the bulk dielectric function, or to its energy 
^ \ derivative, it is growing the belief that surface optical 
O ' spectra are mostly determined by bulk effects, and there- 
fore not very useful as a tool of surface characterization. 
In this paper, we discuss the origin of these bulk-like 
features, and at the same time emphasize the presence 
in optical spectra of other features, more related to the 
C^ ' surface structure. 

In 1996, Rossow, Mantese and Aspnes [p| recognized 
that RAS data on chemically saturated surfaces gener- 
ally resemble the energy derivative lineshapes of the cor- 
responding bulk spectra, dlm[e},{'^)] / d{ijj) , while surfaces 
with unsaturated dangling bonds (DBs) often yield RAS 
lineshapes resembling the bulk spectrum, Im[eb{oj)\. 
They explained the latter lineshapes in terms of surface- 
induced changes of the electron-hole interaction and of 
local fields, while derivative-like spectra were explained 
in terms of surface-perturbations on the energies of bulk 
states. From these findings, they inferred the occurrence 
of shorter lifetimes of electrons and holes near the sur- 
faces. We show here that this deduction is not nec- 
essary. Furthermore, we demonstrate the existence of 
other mechanisms able to produce energy-derivative line- 
shapes. 
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A further step along the way of attributing most RAS 
and SDR features to bulk effects has been done by Uwai 
and Kobayashi (UK) in 1997 Q]. They measured surface 
photoabsorption (SPA) spectra for different conditions 
of the GaAs (001) and (111) surfaces, from which the 
changes of the surface dielectric tensor were extracted. 
The imaginary parts of such changes have peaks at 2.6-3 
eV and at 4.5-4.7 eV, close to the main structures, Ei 
and E2, of the bulk dielectric function. The lineshapes 
are similar to the imaginary part of the bulk dielectric 
function in the case of the (001) surface, and to its deriva- 
tive for the (111) surface. UK conclude that these two 
peaks are not due to transitions involving surface states, 
but to modified bulk electronic transitions. They claim 
that the surface termination effect, first discussed by one 
of the present authors in 1975 M, is responsible for the 
occurrence of bulk-like features in surface spectra. Ac- 
cording to UK, this effect mostly consists in a reduction 
of the polarizability below the surface, arising from the 
quenching of bulk-state wavefunctions near the surface, 
due to their vanishing outside the crystal. This might 
be a long-range effect, extending one hundred Angstroms 
below the surface, which might be hardly included in slab 
calculations. We show here that, although the crystal- 
termination effect is in fact present, the way it has been 
described by UK is rather naive; not a bare reduction of 
the polarizability, but a distortion of its lineshape must 
occur (and indeed occurs), to produce a nonvanishing 
RAS or SDR signal. However, the resulting effect is by 
no means of long range, and is in fact included in slab 
calculations. Moreover, while the crystal- termination 
effect often yields derivative-like lineshapes, we have not 
found bulk-like spectra arising from it. 

The GaAs(llO) surface is a good test case for our calcu- 
lations and discussions, because of its well defined atomic 
structure and for the occurrence of (modest) surface ef- 
fects partially overlapping in energy with (predominant) 
bulk effects [|§. The As-rich GaAs(lOO) /32(2x4) surface 
will also be considered. 



II. THEORY 

We calculate the surface contribution to reflectance, 
that is its relative deviation with respect to Fresnel for- 
mulas, according to the three-layer model 0]: 

Re \ eb[Lu) - 1 / 

for s-light polarized parallel to the i-direction (i = x or 
y) in the surface plane, where 9 is the angle of incidence, 
d the depth of the surface layer, esi(^) is the ii diago- 
nal component of the surface-layer dielectric tensor, and 
eb(ti;) is the isotropic bulk dielectric function. For p- 
light incident in the iz-plane, the anisotropic three-layer 
model yields 0]: 
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The surface-layer dielectric tensor, assumed to be diago- 
nal, is obtained by subtracting the bulk dielectric func- 
tion from the calculated slab dielectric tensor, with a 
suitable choice of the surface-layer depth, d. The re- 
flectivity for s light comes out to be independent of the 
choice of d, coincident with the microscopic formulas not 
relying on the three-layer model. This model is instead 
needed to obtain p-light reflectivity by avoiding the com- 
putationally very demanding inversion of the dielectric 
susceptibility tensor. 

III. RESULTS 

We start by calculating the normal-incidence re- 
flectance anisotropy (RA), 2{Rx — Ry)/{Rx + Ry), of 
GaAs(llO). The latter is the cleavage surface of GaAs, 
and, despite being not reconstructed, undergoes large 
relaxations. Its equilibrium structure, known as the 
"rotation- relaxation model" , is well known both from 
the experimental and the theoretical sides: the sur- 
face As atoms relax toward the vacuum, and Ga atoms 
move in the opposite direction, recovering a quasi-planar 
sp^ bonding with their three As neighbors M. We 
represent the surface using a slab of 31 atomic layers, 
where the actual atomic positions are taken from a Car- 
Parrinello total energy minimization [llO| . Since the slab 
has two equivalent surfaces, the computed slab polariz- 
ability must be divided by two. We then also consider a 
polar surface of GaAs: the As-rich (100) /32(2x4). The 
latter is known to be the stable reconstruction for this 
surface |ll[| , and is a regular array of two As dimers and 
two dimer vacancies (the unit cell contains only two As 
dimers), aligned along the [110] direction. Also in this 
case the actual atomic positions are taken from a Car- 
Parrinello total energy minimization |1J]. In the case of 
GaAs(lOO), since geometry does not allow to build a slab 
with two equivalent surfaces, the calculation is done for 
a system with only one surface reconstructed, i.e. by 
including a real-space cutoff function (a squared cosine, 
approaching one on the interesting surface and zero on 
the other), in the optical transition probability calcula- 
tions, to eliminate the contribution of the back surface 
| p3[ . The electronic states of the slab, as well as of those 
of bulk GaAs, are calculated according to the sp^s* tight- 
binding method, as in Ref. H. The As-As tight-binding 
interaction parameters are those of ref. [Q . The imag- 
inary part of the slab dielectric function is obtained by 
considering transitions at a number of k points in the 
irreducible part of the two- dimensional Brillouin Zone 
(IBZ). The first issue we address is the number of k points 
which are needed to obtain a good convergence. In Fig. 1 
we show the RA of GaAs(llO) calculated with 256, 1024 



and 4096 special k-points in the IBZ; the curves corre- 
sponding to the first two cases are clearly distinct from 
each other. The calculation with 4096 k points, instead, 
is almost coincident with that with 1024 k points. This 
means that 1024 k points are needed to achieve full quan- 
titative convergence of the GaAs(llO) RA. This result 
might be a peculiar property of this and similar surfaces. 
In the case of the GaAs(lOO) /32(2x4) surface, a good 
convergence of the spectrum is already obtained using a 
number of k-points equivalent to 64 in the (1x1) surface 
cell. Similar calculations carried out on Si(110):H show 
that the RA spectrum is already converged with 64 k 
points (usually, calculations are made with 64 k points 
or less [§,^,01, since, even with a well converged k-point 
summation, only qualitative accuracy can be achieved, 
due to the neglect of excitonic and local-field effects.) 

The calculated RAS for GaAs(llO) is qualitatively 
similar to previous calculations, carried out using tight - 
binding or ab-initio methods, and to experiments ||6|,[7|,|l5|. 
The peak at about 2.9 eV embodies a substantial contri- 
bution of transitions across surface states or resonances 
(at variance with Ref. [p|, but in agreement with ref- 
erences [0 and [|5|), while the higher-energy structures 
are essentially due to transitions across surface-perturbed 
bulk states. The main effect of the k-point convergence 
achieved in the present calculation was to reduce the in- 
tensity of the dip just above the 2.9 eV peak and of the 
subsequent structures. 

Having achieved quantitative convergence with respect 
to the number of k points, we can look now at the con- 
vergence with respect to the number of layers, which is 
the main interest here, since changes in lineshapes occur- 
ring for very thick slabs would indicate the presence of 
the long-range effect assumed by UK Q. In Fig. 2a we 
show the GaAS(llO) RA calculated using slabs of 11, 31 
and 93 layers and 1024 k points. The latter two curves 
are almost indistinguishable, while the 11-layer curve is 
also close to them. This means that the calculation has 
already converged with 31 layers, and that the aforemen- 
tioned long-range effect does not occur. The same is 
true for the polar (100) surface: in fig 2c we show the 
calculated RA for the GaAs(lOO) /?2(2x4) surface, where 
small differences show up using slabs of 16, 20 and 40 lay- 
ers. The slow convergence with slab thickness observed 
in calculations W when a smaller set of k points is used 
is therefore due to the error caused by the small number 
of k points, which randomly varies with the number of 
planes. This is nicely demonstrated by Fig. 2b, where 
the same series of slabs as in Fig 2a (11-31-93 layers) 
has been used to compute the GaAs(llO) RAS spectrum 
with a set of 256 k~points. 

In view of the different structures of Eqs. (1), for s 
light, and (2), for p light, one could speculate that the 
long-range effect might cancel in the former case, and 
appear in the latter. To check this possibility, we present 
in Fig. 3 the surface contribution to p-light reflectance 
calculated using the anisotropic three-layer model at an 
angle of incidence of 60 degrees. Again no difference is 



present between the curves calculated with 31 and 93 
layers, definitely showing that surface optical properties 
are well converged with slabs of 31 layers. The present 
results increase our confidence in slab calculation, not 
only since very long-range effects, which can be hardly 
embodied therein, are excluded, but also because thin 
slabs, as the 11 and 16-layers ones, which are the only 
ones that can be afforded in ab-initio calculations 0, 
already yield rather good results. 

Let us discuss now in more detail the crystal- 
termination effect invoked by UK 0] . It has been firstly 
addressed by one of the present authors in connection 
with the reflectivity at the direct minimum gap of a semi- 
conductor g . By disregarding the microscopic structure 
of the surface, and describing it just as an infinite po- 
tential barrier preventing electrons from escaping into 
vacuum (the crystal termination) , the reflectivity was ob- 
tained starting from the wavefunctions calculated accord- 
ing to the effective-mass approximation. Since the wave- 
functions must vanish at the crystal-termination plane, 
the envelope plane-waves occurring in an inflnite crystal 
are replaced by sine-type standing waves. When look- 
ing at the imaginary part of the local dielectric function, 
Im[e{z,u))], this yields a region below the surface where 
this quantity is smaller than in the bulk crystal. The 
depth of such a region is of the order of ir/kz , k^ being the 
largest wavevector of the relevant transitions. In the case 
of the minimum gap, hence, fc^ — [2m*{hLj — Eg)/h'^]^/^, 
where m* is the reduced electron-hole effective mass and 
Eg the direct- gap energy. By taking an effective mass 
of 0.1 and hui — Eg as 0.1 eV, we estimate this depth to 
be of the order of 60 Angstrom. This is the quenching of 
the (bulk) dielectric function below the surface, that UK 
assume as the most important effect. However, crystal 
termination affects the optical properties also in another 
way: since the matrix elements of the momentum op- 
erator must be calculated between the surface-perturbed 
wavefunctions (sine- type in the case discussed here), they 
are different from those of the infinite crystal. More 
explicitly, kz is no longer conserved in subsurface opti- 
cal transitions; the breaking of this selection rule yields 
spectra distorted with respect to (namely, broader than) 
the corresponding bulk spectra, acting as an additional 
broadening localized near the surface. We will show be- 
low that this additional broadening is the most important 
crystal-termination effect infiuencing the surface optical 
properties of GaAs. 

Differently from the case discussed above, the main 
structures of bulk spectra, which yield the most promi- 
nent bulk-related structures in RAS and SDR, are due to 
transitions at saddle-points of the joint density of states. 
The characteristic k^ at saddle points is much larger than 
at the direct gap, because a large region of Brillouin Zone 
is available for optical transitions at the saddle-point en- 
ergy. For instance, vertical transitions along all the A line 
are responsible for the Ei structure in GaAs. Hence the 
largest k^ is of the order of the BZ boundary, n/a, and 
the depth of the surface-perturbed region is of the order 



of the lattice constant, a. This explains why we have 
not found in Figs. 2 and 3 any indication of long-range 
effects close to saddle- point energies. 

The simplest model of quenching is to assume that the 
polarizability is completely suppressed within some depth 
d below the surface. This, however, would be equivalent 
to shift the surface by d, and would not give any contribu- 
tion to the reflectance. Hence we consider a slightly dif- 
ferent model, where the polarizability is partly quenched, 
say 50 percent, in a depth d. In practice we assume, 
within the depth d, a surface dielectric function of the 
form: 

esi{uj) = ii ■ Sb{Lo - ALUi,"/i) (3) 

where £6(w) is the bulk dielectric function, f^ (< 1) repre- 
sents the quenching, ji is the broadening (possibly differ- 
ent from the bulk one), and Auji is a possible frequency 
shift. When Acji = 0, 7^ = ^buik, and U = 1, esi(w) 
coincides with the bulk dielectric function. Taking fi < 1 
with 7i = Jbuik and Auji ~ would not modify the s- 
light reflectivity, since the numerator and denominator 
in equation (1) are proportional to each other, and hence 
the fraction is a real number, with vanishing imaginary 
part. However, this is not the case for p-light reflectiv- 
ity, which may undergo some change. The full line in 
Fig. 4a shows the surface contribution calculated in this 
way. It is clear from the figure that this model has no 
relation with the output of the slab calculation (dashed 
line); hence the pure quenching effect cannot account for 
the surface contribution to reflectance. 

We consider next the pure broadening model, i.e. 
fj = 1, AiOi = 0, and 7^ > jtuik- Now the surface is 
assumed to have the same dielectric function as the bulk 
has, but with broader lineshapes, as a consequence of 
the breaking of the fcx-conservation near the surface. In 
Fig. 4b we show the surface contribution to the reflec- 
tivity of normally incident light, with polarization per- 
pendicular and parallel to the [110] chains, as calculated 
from the slab polarizability (dashed and dotted lines), 
and according to the broadening model (full line). We 
assume a broadening of 100 meV at the surface, while 
it is 30 meV in the bulk. The curves are rather simi- 
lar, showing lineshapes resembling the energy-derivative 
of the imaginary part of the dielectric function. Also the 
surface-state related peak at about 2.9 eV is embodied 
in the broadening-model spectrum (this occurs only be- 
cause the peak mentioned above overlaps in energy with 
the El bulk structure around 3 eV). However, the dif- 
ferences between the broadening model and microscopic 
calculations, which seem to be small in this spectrum, 
become very large in the RA spectrum, shown in Fig. 
4c. Here the dashed line is obtained as in Fig. 1, that 
is from the 31-layer slab calculation. We can produce a 
RA-curve according to the broadening model by assum- 
ing that the depth where the dielectric function is broader 
than in bulk is different for the two polarizations, or, in 
an equivalent manner, that also some quenching of the di- 
electric function occurs (f^ < 1), whose amount depends 



on the direction of light polarization {i^ 7^ iy)- By assum- 
ing a suitable depth or quenching difference, we produce 
the full line in Fig. 4c, which is of course proportional to 
the continuous line in Fig. 4b. The two curves in Fig. 4c 
are markedly different, although the peak at 2.9 eV (the 
only spectral feature related to transitions across surface 
states!) is present in both curves. 

Hence we can conclude by recognizing the occurrence 
of bulk- derivative- like features in surface optical spectra 
calculated for a given polarization of light, due to the 
broader lineshape of the dielectric function near the sur- 
face. This broader lineshape is due to the breaking of the 
kz conservation rule (namely, it is a crystal-termination 
effect), is included in slab calculations, and does not im- 
ply a shorter lifetime of electrons near the surface than 
in bulk. When anisotropy difference spectra are taken 
for GaAs(llO), however, these features largely cancel, so 
that the surviving RA has about no relation to the broad- 
ening model. Of course, such cancelation may be smaller 
at other surfaces, so that derivative-like lineshapes may 
be present in RAS and SDR spectra. 

As a last point, we can assume that the (bulk) dielec- 
tric function near the surface can undergo small shifts 
of peak positions {Auji ^ 0), in addition to broadening 
and quenching. To this aim it is not needed, as assumed 
in Ref. 1^, that electrons and holes excited in optical 
transitions are kept close to the surface by their short 
life-times, in order to be shifted in energy by the sur- 
face potential. The required small shifts of the peaks 
of the surface dielectric function may be produced by 
the surface-perturbation on the wavefunctions and, con- 
sequently, on the local polarizability. It is a matter of fact 
that the layer-projected density-of-states may be differ- 
ent from the bulk one. The same, of course, can occur 
for the z-dependent dielectric function [|l6|, whose aver- 
age over the first few layers yields the surface dielectric 
function. 

We have tried to obtain bulk-like difference (RAS or 
SDR) spectra by suitably varying f^, 7^, and Awi, i.e. 
by shifting, broadening and quenching the bulk spec- 
trum. By varying the parameters above, we often ob- 
tained derivative-like spectra, never obtained bulk- like 
spectra, and sometimes hybrid spectra (see Fig. 5a, full 
line). It is worth to notice that this hybrid spectrum is 
rather similar (although energy shifted) to the microscop- 
ically calculated RAS spectrum, also shown in Fig. 5a 
(dashed line). For some choice of the parameters we got 
difference spectra approximately bulk like, that is showing 
peaks close to the two bulk critical-point energies, but, 
differently from the bulk spectrum, with a negative region 
in between (Fig. 5b, full line). A similar RA spectrum 
is the result of a realistic TB slab calculation carried out 
for another GaAs surface, the (polar) (100) /32(2x4). The 
calculated RAS is shown in Fig. 5b, by the short-dashed 
line, while the experimental spectrum, more similar to 
the bulk one, corresponds to the long-dashed line p7[ . 
This suggests that the surface-exciton and surface local- 
field effect may be determinant to yield bulk-like surface 



spectra. A recent calculation for Si(110):H [18|, where 
the experimental RA lineshape is bulk like [19|, shows 
indeed that the surface local- field effect, treated therein 
according to the point-dipole approximation, is crucial to 
obtain a bulk-like theoretical lineshape. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we have shown that surface effects on 
optical properties of GaAs are localized in a few mono- 
layers below the surface. As it has been discussed in 
section III, these results do not depend explicitly on the 
particular system considered, suggesting a more general 
validity, i.e. indicating that for a wide class of semi- 
conductor surfaces the surface effects on optical prop- 
erties, including the crystal termination effect, are well 
described using slabs of a few tens of monolayers. In 
the absence of peculiar features due to surface states, 
the crystal- termination effect can be phenomenologi- 
cally modeled as a shift, broadening and reduction of 
the bulk spectrum. Many combinations of these param- 
eters yield surface spectra resembling the derivative of 
the bulk absorption spectrum. It must be emphasized 
that the amounts of shift, broadening and reduction are 
ultimately determined by the microscopic structure of 
the surface; furthermore, these bulk-derived structures 
coexist with transitions directly involving surface states. 
After subtraction of individual spectra to obtain RAS or 
SDR spectra, the resulting lineshape can be qualitatively 
different from a derivative- like lineshape, as in the case of 
GaAs(llO). On the other hand, approximately bulk-like 
spectra are obtained for some values of the parameters. 
However, truly bulk-like spectra can hardly be obtained 
in terms of the crystal-termination effect, and they did 
not even occur as results of our realistic slab calculations. 
Hence, many-body effects like the surface-exciton or the 
surface-local field effect, not included in the one-electron 
theory, seem to be determinant to obtain truly bulk-like 
lineshapes. 

In conclusion, we agree with Rossow et al. |^ and with 
UK Q that some features of surface spectra originate 
from transitions across bulk states. Surface termination 
effects, however, involve a more complex mechanism than 
that described by UK. In fact, the broadening of the bulk 
dielectric function near the surface is the most important 
crystal termination effect, due to the breaking of kz con- 
servation at surfaces. It does not imply, however, that 
photogenerated electrons and holes have shorter lifetimes 
than in the bulk. This effect yields derivative-like line- 
shapes, as those obtained at many chemically-saturated 
surfaces. 

Finally, we would like to stress that many effects con- 
cur to determine surface optical properties. It is not pos- 
sible to interpret optical spectra of all surfaces in terms 
of a single effect, either the crystal termination effect, or 
transitions across surface states. Caution must also be 



used in assigning spectral features to bulk-state transi- 
tions uniquely because of their energy positions, as ex- 
emplified by the case of GaAs(llO), where we found that 
the main peak of the calculated spectrum, occurring al- 
most at the same energy as the Ei bulk feature, contains 
a substantial contribution of transitions across surface 
states. 
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FIG. 1. Calculated Reflection Anisotropy Spectrum of 
the GaAs(llO) surface, and its convergence with respect to 
the Brillouin zone sampling. Full line: 4096 k-points in the 
irreducible wedge of the Surface Brillouin zone; long-dashed 
line: 1024 k-points; dotted line: 256 k-points. In the latter 
case, full convergence has not yet been reached. 



FIG. 2. a) Convergence of the theoretical RA spectrum 
of GaAs (110), calculated using 1024 k points, with respect 
to the thickness of the slab (number of atomic layers) used 
in the calculation. Full line: 93 layers; long-dashed line: 31 
layers; short-dashed line: 11 layers. The results for 31 and 91 
layers are almost identical, b) The same, calculated using 256 
k points, c) Convergence of the theoretical RA spectrum of 
GaAs (100) /32(2x4), calculated using 8 k points (equivalent 
to 64 in the (1x1) surface cell), with respect to the thickness 
of the slab. Full line: 40 layers; long-dashed line: 20 layers; 
short-dashed line: 16 layers. 



FIG. 3. Surface contribution to reflectance for p-light 
incident in the y-z plane at 60 degrees, for GaAs(llO) slabs 
of 31 layers (dashed line) and 91 layers (full line), using 1024 
k-points in the SBZ. The y direction is parallel to the chains 
in the surface plane, while z is perpendicular to the surface. 
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FIG. 4. a) Microscopic calculation (110 slab) of the sur- 
face contribution to reflectance for p-light incident in the y-z 
plane at 60 degrees (dashed line), compared with the results 
of the "quenching model" (full line); b) Microscopic (slab) cal- 
culation of the surface contribution to reflectance for normally 
incident light polarized along the (110) chains, ( dotted line) 
and perpendicularly to them ( dashed line), computed with 
a 31-layers GaAs(llO) slab and 1024 k-points in the ISBZ, 
in comparison with the results of the 'broadening model' (full 
line); c) "Broadening model" results for the RA spectrum of 
GaAs (full line), compared with the microscopic calculation 
of Fig. 1 



FIG. 5. a) full line: GaAs RAS results from the quench- 
ing-broadening-shifting model, with the following choice of 
parameters: i^ = 1, f^ = 0.5, jx = 0.1 eV, jy = 'jtuik = 
0.03 eV, Aujx = Auiy = O.leV (see text). Dashed line: mi- 
croscopic calculation for the (110) slab, with 31 layers and 
1024 k-points. b) long-dashed line: experimental RAS data 
for GaAs (100) /32(2x4), from ref. 0. Short-dashed line: mi- 
croscopic calculation for the same surface. Full line: quench- 
ing-broadening-shifting model, with the following choice of 
parameters: ix = iy = 1, 7x = 7y = O.leV, jbuik = O.OSeV, 
AiUx = -O.leV^, Aojy = 0. 
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